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Abstract. We consider the issue of non-Mar kovianity of a quantum dynamics 
starting from a comparison with the classical definition of Markovian process. We 
point to the fact that two sufficient but not necessary signatures of non-Markovianity 
of a classical process find their natural quantum counterpart in recently introduced 
measures of quantum non-Markovianity. This behavior is analyzed in detail for 
quantum dynamics which can be built taking as input a class of classical processes. 
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1. Introduction 

In the field of open quantum systems the interaction with an external environment 
introduces a stochasticity element in the dynamics, which has typically been described 
as a quantum process in analogy with classical processes [1]. Quantum Markovian 
processes have been understood as described by a master equation in Lindblad form 
[2, 3], which grants complete positivity. Furthermore these dynamics can be understood 
as an average over trajectories corresponding to suitable measurements continuous in 
time on the side of the environment, thus providing a direct link with the description 
of quantum measurement [4]. 

However the very notion of non-Markovianity for a quantum dynamics has long 
remained vague, being usually associated with the occurrence of revivals or non 
exponential relaxation, apart from mathematical work in which a proper definition 
of a quantum process has been addressed [5, 6]. Most recently, different approaches 
have been devised in order to assess and quantify non-Markovianity in the dynamics 
of an open quantum system, looking at states rather than at multitime correlation 
functions. One approach, based on the idea that memory can be seen as a backflow of 
information from the environment to the system, studies the distinguishability among 
states in the course of the dynamics, identifying non-Markovianity with revivals in the 
distinguishability. The amount of non-Markovianity then depends on frequency and 
relevance of these revivals [7]. Note that this approach is actually valid in order to 
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estimate non-Markovianity of the dynamics in any statistical theory, simply considering 
the mathematical representation of the set of states and a distance on it, which should 
be a contraction under the action of positive probability preserving maps. Another 
approach relies on the divisibility over arbitrary intermediate time intervals of the 
evolution map into well defined maps preserving positivity and probability. It estimates 
non-Markovianity quantifying the failure of this property [8] . 

In recent work it has been pointed out that these two approaches can be connected 
to certain signatures of non-Markovianity to be detected at the level of the one-point 
probability in a stochastic process [9], and a few examples of this behavior have been 
considered building on a class of stochastic processes known as semi-Markov processes 
[10, 11]. In this contribution we will consider the class of maps which can be built 
within this framework by considering quantum dynamics determined by a generic 
unital stochastic map and an arbitrary waiting time distribution. We will exploit these 
examples to study in detail the different performance of the two measures and the variety 
of possible behavior which can be obtained. This class of dynamics will also allow us 
to point to the proper relationship between divisibility of the time evolution in terms 
of completely positive maps and behavior of the coefficients in the time-convolutionless 
form of the equations of motion. 

The paper is organized as follows. In Sect. 2 we introduce the definition of classical 
Markov process, together with previous proposals to extend this definition to the 
quantum realm. We further consider two sufficient signatures of non-Markovianity to 
be read from the one-point probability density. In Sect. 3 it is shown how the extension 
of these signatures to the quantum case recovers two recent proposals for the definition 
of a non-Markovian dynamics, which are studied and compared for a general class of 
examples. Conclusions and final remarks are presented in Sect. 4. 

2. Markovian and non-Markovian processes 

The notion of stochastic process is of the greatest importance in the description of 
phenomena in which randomicity appears. Due to interaction with the environment 
open quantum systems do provide a natural setting in which a stochastic description 
appears. This in turns leads to the difficult question about how to properly characterize 
a quantum stochastic process, addressed in different ways by physicists [12, 1] and 
mathematicians [5, 6]. Within this framework a further important question is how 
to define Markovianity of a process and how to associate a memory to the various 
phenomena. 

In the classical case one has a precise definition of Markovian process, relying on the 
knowledge of all finite dimensional distributions of the process, which is best formulated 
in terms of conditional probabilities [13] 

Pl\n{pm tn\%n—li tn—li ■ ■ ■ i ^o) Pl\\{%m tn\%n— It tn—l) (-Q 

with t n ^ t n -i ^ . . . ^ ti ^ to. It tells that the probability for the random variable 
to assume the value x n at time t n only depends on the last assumed value, and not 
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on previous ones, thus properly formalizing the notion of lack of memory. The set of 
conditions Eq. (1) lead in particular to the Chapman-Kolmogorov equation obeyed by 
the two point conditional probability, also called propagator 

p 1 \ 1 (x,t\y,s)= ^pi|i(x,^,r)pi|i(z,r|y,s), (2) 

z 

valid for t ^ r ^ s. A solution of this equation fully characterizes a Markovian process 
in that all finite dimensional distributions can be obtained given the initial probability 
distribution. Since the propagator is on its turn determined from the dynamics of all 
mean values, this means that the whole process is known from the dynamics of mean 
values, a result known as regression theorem [14]. All multitime correlation functions 
of the process are thus fixed from the propagator and therefore from the mean values. 
Indeed Lindblad introduced a definition of quantum Markovian process by relying on 
the validity of the quantum regression theorem [5]. Consider a system S interacting 
with an environment E according to a unitary evolution U (t) and multitime correlation 
functions of the form 

(M n (t n ) . . . Mi(h)) = Tr s Tr E {M n U (t n - t^Mn-i . . . MiU^ps ® Pe}, 

where M.i denotes a quantum operation corresponding to measurement of a certain 
system observable, e.g. A4(g) = Ylk^ePk for the von Neumann measurement of the 
self-adjoint operator A = Y2k a k-Pk- If the quantum regression formula applies such 
multitime correlation functions can be expressed as 

(M n (t n ) . . . Mi(ii)) = Tr S {M n $(t n - tn-^Mn-! . . . Mi^(h)p S }, 

with $(£) completely positive maps acting on the system only. In the weak and singular 
coupling limit these maps can be shown to satisfying a semigroup composition law 
[15], so that they can be expressed with a generator in Lindblad form. As a result in 
particular mean values and correlations do obey the same dynamical equations. This 
result, though by no means obvious for the behavior of the reduced dynamics of an open 
system, is a basic working tool in open quantum system theory [1], and the extension of 
the validity of the regression formula to the non-Markovian setting is an issue of great 
relevance. The problem has been recently addressed for a special class of non-Markovian 
evolutions known as generalized Lindblad structure [16, 17, 18], showing however that 
regression formula are generally not valid [19]. 

2.1. Witnesses and quantifiers of non-Markovianity 

The validity of the Chapman-Kolmogorov equation for a Markovian process entails 
two simple consequences on the behavior of the propagator, which are at the heart of 
two recently introduced quantum measures of non-Markovianity [7, 8]. Let us consider 
for the sake of simplicity a random variable taking values on a finite set, so that the 
propagator from time s to time t can be written as a stochastic matrix A(t, s) and the 
state of the system at a given time t can be expressed as a probability vector pit). As 
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distance among probability vectors we consider the Kolmogorov distance [20] 
D K (p\t),p\t))= 1 -J2\pK t )-pM- 

i 

Upon validity of the Chapman-Kolmogorov equation one has p(t) = A(t, s)p(s) for 
any ^ s ^ t, and exploiting the definition of stochastic matrix this entails that 
the Kolmogorov distance among two probability distributions following a Markovian 
dynamics do decrease monotonously in time 

D K (p\t),p 2 (t)) < D K (p 1 (s),p 2 (s)) Vt > 8 . 

This relation naturally provides a necessary condition for the Markovianity of a 
stochastic process, stating that in a Markovian process the one-time probability densities 
or distribution functions evolving from two distinct initial situations do get less and less 
distinguishable. A non monotonic behavior in time of the Kolmogorov distance among 
two states thus provides a witness of the non-Markovianity of the process. Of course 
this condition now appears as a sufficient condition only, or if one prefers a different 
definition of Markovianity. 

Let us stress that the Kolmogorov distance actually corresponds to the L\ distance 
among probability vectors, thus being a natural distance in any statistical theory. Indeed 
a statistical theory relies on the existence of two spaces, one dual to the other, in which 
states and observables do live. The specific choice of spaces and their commutativity 
or non commutativity do fix the statistical structure of the theory [21, 22, 23, 24]. In 
the case of classical mechanics one considers the L\ space of probability distributions 
on phase space, while the observables are given by the space of bounded functions. 
Considering quantum mechanics the dual Banach spaces are given by the space of trace- 
class operators with the trace norm topology, in which states are described by statistical 
operators, and the space of bounded operators with the uniform norm, in which to 
consider the observables of the theory. The Kolmogorov distance thus naturally becomes 
the trace distance among statistical operators. 

The validity of the Chapman-Kolmogorov equation for a Markov process entails 
another important consequence at the level of the one-point probability density, namely 
rewriting Eq. (2) in terms of stochastic matrices for a finite dimensional process one has 
A(t, s) = A(t, t)A(t, s) for any s ^ r ^ t. Taking s as an initial time set equal to zero 
one has the relation 

A(t,0)=A(t,r)A(r,0) (3) 

where the crucial fact is that each A is a well-defined stochastic matrix sending any 
probability vector to a probability vector. For a Markov process these matrices are 
fixed by the transition probabilities. Eq. (3) then describes a divisibility property which 
provides a witness of Markovianity. However the violation of this property only provides 
a necessary but not sufficient condition in order to ascertain Markovianity. Indeed also 
for non-Markovian processes one might find a collection of stochastic matrices obeying 
Eq. (3), which nevertheless do not coincide with the transition probabilities of the 
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process [25]. Otherwise stated these collection of stochastic matrices can be taken as 
the transition probabilities of a Markov process which at the level of the one-point 
probabilities cannot be distinguished from the original one. 

These signatures of non-Markovianity for a classical process correspond to criteria 
used to define a non-Markovian dynamics in the quantum case in [7] and [8], by 
considering distinguishability as quantified by the trace distance, and divisibility as 
characterized by a composition in terms of completely positive maps respectively. In 
order to further assign a weight to the deviation from the Markovian behavior, relying 
on these signatures two measures of non-Markovianity have been introduced, which 
essentially assign a weight to the time intervals in which either trace distance grows or 
divisibility, at the quantum level to be understood as divisibility in terms of completely 
positive maps, fails. These weights can actually be assigned in different ways. For 
example one can consider quantifiers of distinguishability on the state space other than 
trace distance, provided they are contractions with respect to the action of positive 
maps [26, 27], e.g. relative entropy which appears in a natural way in certain dissipative 
systems [28, 29]. 

2.2. Semi-Markov processes 

Let us now consider a class of stochastic processes for which a simple characterization 
is available, and which include both Markovian and non-Markovian processes, namely 
semi-Markov processes [30] . Such processes arise considering transitions among different 
sites determined from certain jump probabilities, fixed by a stochastic matrix as in 
a Markovian chain, and random waiting times between jumps determined by site 
dependent waiting time distributions. They are characterized via a semi-Markov matrix 



whose entries give the probability density for a jump from site n to site m in a given time 
t. The transition probabilities 7r mn build up a stochastic matrix, while the f n (j) are the 
waiting time distributions which actually characterize whether the process is Markovian 
or not according to the classical definition. Indeed such a process is Markovian only if 
all waiting time distributions are memoryless, that is exponentially distributed 



In all other cases such processes are non-Markovian. This simple characterization of 
semi-Markov processes, together with the fact that their transition probabilities can be 
obtained as solution of a generalized master equation with a memory kernel, allows to 
exemplify the meaning of the signatures of non-Markovianity introduced above. In the 
case in which all waiting time distributions are equal, though otherwise arbitrary, the 
transition probability T(t, 0) of such a process obeys a closed integrodifferential equation 
given by 



(Q)mn(r) = (n) mn / n (r), 



(4) 



fn(r) = \ n e 



(5) 




(6) 
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where k(r) is a memory kernel fixed by /(r) according to the relationship 

/(r) = f dtk(r - t)g(t), g(t) = 1 - f dr/(r), (7) 

where g(t) is the survival probability associated to /(t). 

We now consider a two-dimensional system and take two equal waiting time 
distributions, so that the semi-Markov matrix of Eq. (4) can be written 

with 7r, a - the jump probabilities from one site to the other. Starting from Eq. (6) for 
the case n — a — 1/2, one can express the transition probability as a function of the 
survival probability only, and using T(t,s) = T(t, 0)T~ 1 (s, 0) one obtains the matrices 

T( , v i ( i + g(t)/g(s) i - <K*)M«) \ (q) 

[h ' 2 ^ 1 - <Kf)/</(a) l + g(t)/g(s) y W 

which connect probability vectors at a time s with probability vectors at a later time t. 
Since g is a survival probability, these matrices are stochastic matrices for any couple 
t ^ s, independently on the fact that the associated semi-Markov process is Markovian 
only if / is exponentially distributed. Thus for any choice of /(r) different from the 
exponential one has an example of process which is non-Markovian in the classical sense, 
but still whose one-point probabilities do contract in a monotonous way with respect to 
the Kolmogorov distance and are connected by stochastic matrices. This indeed shows 
that contractivity under the Kolmogorov distance and divisibility are sufficient but not 
necessary criteria in order to detect non-Markovianity in the classical sense. 

As a complementary situation, let us consider the case ir — 0, a — 1, so that once 
in a state the system jumps with certainty to the other, thus obtaining 

q(t)/q(s) l-q(t)/q(s) 
q(t)/q(s) l + q(t)/q(s) 




(10) 



The role of the survival probability g(t) is here replaced by the quantity q(t) whose 
Laplace transform reads 

1 1 - f(u) 
qiu) = 7. . 

«! + /(«) 

Recalling that the probability for n jumps in a time t for a waiting time distribution 
fit) is given by 

p n (t)= [ drf(t-r)p n ^(r), (11) 
Jo 

so that p n (u) = po(u)f n (u), one finally has 

oo oo 

q(t) = J>2n(t) " J>2n+l(0- (12) 

n=0 n=0 
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Figure 1. (Color online) (top left) Kolmogorov distance for a classical semi-Markov 
process characterized by ir = and a — 1, and a waiting time distribution given by 
the convolution of two exponential ones with A2 = 0.5Ai. The different trajectories 
correspond to initial probability vectors such that |p} (0) — (0) | are equally distributed 
between and 1. Revivals clearly appear, (bottom left) Kolmogorov distance for a 
semi-Markov process with the same waiting time distribution but tt = a = 1/2. The 
distance decreases monotonically. In both cases the process is non Markovian according 
to the classical definition, (right) The modulus of the functions q m , corresponding to 
the convolution of m equal exponential distributions, for m up to 6 and in the case of 
a semi-Markov process with 7r = and a = 1. The horizontal lines denote the height 
of the first non trivial maximum of \q m (t) |, quantifying the non-Markovianity of the 
first interval in which trace distance grows. 

The quantity q(t) therefore expresses the difference between the probability to have 
an even or an odd number of jumps. It immediately appears that this quantity is 
not necessarily positive, so that the matrices defined by Eq. (10) are not necessarily 
stochastic matrices, and apart from the case of an exponential waiting time distribution, 
corresponding to a truly Markovian process, different expression of f(r) may or may 
not lead to contractivity under the Kolmogorov distance and divisibility. These 
possible behaviors are considered in Fig. 1, where the Kolmogorov distance for two 
non-Markovian classical process is considered. 

3. Quantum non-Markovian dynamics 

The characterization of well-defined classes of quantum time evolutions which have the 
property of being completely positive, though not in Lindblad form, is a highly non 
trivial task. A whole class of completely positive dynamics can be obtained considering 
a quantum generalization of the classical semi-Markov processes introduced in Sect. 2.2, 
also allowing for the connection and comparison with a classical process. Such dynamics 
are given by the solution of integrodifferential equations with a memory kernel which 
is formally of Lindblad type [10, 11]. For the case in which the dynamics allows for a 



8 



clearcut probabilistic reading such equations can be written in the form [16] 

j t p(t)= j\rk{t-r)[£-l]p{r). (13) 

Indeed if the function k(t) can be associated to a well defined waiting time distribution 
f(t) according to the relation Eq. (7), then the time evolution of the solution p(t) of 
Eq. (13) can be expressed as repeated actions of the completely positive and trace 
preserving map £, also called stochastic map, distributed in time according to the 
renewal process fixed by f(t) 

oo 

A(f,0)p(0)= X>n(*)£>(0), (14) 

n=0 

where p n (t) is defined in Eq. (11). Note that Markovianity or non-Markovianity of 
this quantum time evolution, according to either criterion of Sect. 2.1, at variance 
with the classical definition depend on both elements of the couple {£,/(£)}• More 
precisely, while the only truly non-Markovian waiting time distribution corresponding 
to the exponential leads to a delta correlated kernel in Eq. (13), and therefore to a 
Lindblad equation, a generic waiting time distribution might still lead to a Markovian 
dynamics depending on the completely positive trace preserving map £ as we shall see 
below. 



3.1. Unital stochastic maps 

Let us first leave / (t) unspecified and consider as £ the socalled Pauli channel 

£ [P] = \ P + \ X 0- X P<J X + XyO-yPCTy + \ Z (T Z p(T Z , 

with A = (A , \ x , \, X z ) a probability distribution, \ ^ and ^\ A, = 1. As discussed 
in [31] up to unitary transformations this expression provides the most general unital 
stochastic map, which preserves both trace and identity for H = C 2 . The map A(t, 0) 
associated to £ through Eq. (14) can be expressed in terms of a basis of superoperators 
acting on operators in C 2 . The standard representation is 

3 

A[p] = ^kir k Trs[r}pi (15) 

fc,/=0 

with {T fc } fe=0jli2) 3 a basis of operators on C 2 and A M = Tr s {rlA[Ti]}. The most 
convenient choice of orthonormal basis corresponds to Tj = ^gCj, with ao — 1 and 
<7j the standard Pauli matrices. These basis elements are eigenoperator of the Pauli 
channel. The eigenvalues are given by fi = AX, with A 0i = A i0 = 1 and A^ = 25j k — 1 
for j, k — 1, 2, 3, so that p,o — 1 and — 1 ^ pi ^ 1. In order to determine the map A (t, 0) 
we consider its action on the basis elements, given by 

A (t, 0) [<7i] = Xi (t) a*, (16) 

where the quantities Aj (t) are defined as 

oo 

\(t) = X>«(*K> ( 17 ) 

n=0 
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and thus correspond to the generating function of the discrete probability distribution 
p n (t) evaluated at pi, in particular Ao (t) = 1, while — 1 ^ Aj (t) ^ 1. Note that according 
to Eq. (16) the standard representation of statistical operators on the Bloch sphere the 
action of the map A (t, 0) transforms the surface of the sphere into ellipsoids whose axes 
have a time dependent length |Aj (t)\. In the classical setting Markovianity of the semi- 
Markov process is only obtained for an exponential waiting time distribution, so that 
the p n (t) are a Poisson distribution. In this case the generating function reads 

\(t) = e -^ xt , (18) 

where A is the parameter of the exponential distribution, thus corresponding to ellipsoids 
whose axes shrink monotonously according to an exponential law. Exploiting Eq. (11) 
the Laplace transform of the A« (t) can be directly expressed in terms of the Laplace 
transform of the waiting time distribution according to 

= ^4rv < 19) 

By linearity Eq. (16) determines the map, which is fixed by the matrix elements 
Pn(t) =^[l + A z (t)( Pll (0)-poo(0))] 

Pio (t) = ^ [Pio (0) (A, (t) - \ y (*)) + poi (0) (A, (t) + \ y (*))] . 

This provides the general solution of Eq. (13) for arbitrary / (t) and £ a generic unital 
stochastic map. For later use it is also convenient to express the map in a different 
superoperator basis 

3 

A[p] = KiTkwrl (20) 

k,l=0 

with A' kl = J2i^ r s{ T i r j T kM T i\} ■ This representation for superoperators was introduced 
in [32], and apart from a multiplicative factor corresponding to the space dimension it 
associates to the map its Choi matrix. For the previous choice of basis set the associated 
matrix is still diagonal, so that 

A (t, 0) [p] = p (t) p + p x (t) (T x pa x + p y (t) Oypay + fj, z (t) a z pa z , (21) 

with pL {t) = A~ 1 \ (t) = \A\ (t) and the {pi (£)} at any time a probability distribution. 

3.2. Analysis of non- Markovianity 

We now study Markovianity of these quantum maps according to the criteria introduced 
in [7, 8]. The first criterion is based on the behavior of the trace distance among 
distinct initial states, which quantifies how their distinguishability varies in time. For 
the considered map the trace distance reads 

D(p\t),p*(t)) = 1 -\\p\t)-p 2 (t)\\ 1 



XI (t) A p (0) 2 + XI (t) KA c (0) 2 + XI (t) 3A C (0) 2 , 
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where A p (0) and A c (0) denote population and coherence differences at the initial time. 
The time derivative of this quantity, which detects non-Markovianity identified with 
growth of trace distance in time for at least a couple of possible initial states, is 

1>2 \ z (t) X z (t) A p (0) 2 + \ x (t) A, (t) ^A c (0) 2 + A, (t) \ y (t) 5A C (0) 2 

it) A p (0) 2 + A 2 (t) MA c (0) 2 + A 2 (t) 3A C (0) 2 

so that a necessary and sufficient condition for non-Markovianity is that at least one of 
the functions |Aj (t)\ grows in a certain time interval. This condition coincides with the 
requirement of divisibility of the quantum map in terms of positive (but not necessarily 
completely positive) maps according to Eq. (3). For this time evolution the map A (t, s) 
is represented according to Eq. (15) by the matrix 

* U s A - Mf) M*) A,(t) 

A kl (t, s) = dtag 1, — ^, — ^ 
V \ x (s) X y (s) X z (s) / 

which corresponds to a positive map if and only if all eigenvalues lie between zero and 

one. Let us now spell out these general relations for specific choices of {£, fit)}- 

For A = (0, 0, 0, 1) we recover the phase-flip channel 

which does not have a classical counterpart. For this case one has n = (1,-1,-1,1), 
leading to A (t) = (1, q (t) ,q (t) , 1), with q (t) as in Eq. (12). The map A (t, 0) describes 
pure dephasing, with coherences multiplied by a factor q(t). As a result the time 
derivative of the quantifier of the distinguishability among the two evolved states grows 
whenever \q(t)\ grows in time. In such time intervals, which we denote collectively by 
^+ = Ui( a i,bi), we have 

^^(O))^ 1^(0)1^19(01 

so that the maximal growth is obtained for opposite equatorial states on the Bloch 
sphere. The measure of non-Markovianity based on trace distance is then given by 



m)= [ ^^1^)1 = ^(19(^)1-19(^)1) 



To consider an interesting class of situations we now specify also the waiting time 
distribution, considering the convolution of m exponential waiting time distributions 
leading to so called Erlang distributions of the form f m (t) = A ^-i); e ~ Xt , fc> r which the 
difference of the probabilities to have an even and an odd number of jumps becomes 

oo m— 1 



q m (t)=e-*Y,(- 



(Xt) mn+k 



' (mn + k)\ 

For the memoryless case m = 1 one has the strictly monotone decreasing function 
qi(t) = exp (—2Xt) according to Eq. (18), while for m = 2 one obtains 

q 2 {t) = e" A '[cos(At) + sin(At)], 

which oscillates and crosses zero at isolated points. For m ^ 2 these functions exhibit 
an oscillating behavior, so that the minima of |? m (£)| lie on the real axis. The modulus 
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of these functions for values of m up to six is shown in Fig. 1. For each choice of waiting 
time distribution f m (t) the measure of non-Markovianity is given by the series 

A/-(A m ) = \9miU)\, (22) 

UeM 

where M is the denumerable set of times corresponding to the maxima of the function 
\q m \. An exact analytic evaluation of the measure is feasible for m — 2, since the maxima 
correspond to At = nn, with n a positive integer, so that 

1 



A/-(A 2 ) = E e " 



e n - 1 



n=l 

More generally as shown in Fig. 1 the first maximum of \q m \ is above the first 
maximum of \q n \ for m > n, and the same occurs for the other maxima, whose values 
decrease exponentially. This substantiates the statement that waiting time distributions 
corresponding to the convolution of a higher number of exponential terms have stronger 
memory [33]. Indeed in such a case the overall waiting time is the sum of a number 
m of waiting times, which can be thought in series. While each of them is described 
by a memoryless probability distribution, the overall distribution deviates from the 
memoryless case, the more so the higher the number of terms. Another situation 
in which Markovianity or non-Markovianity can be observed is given considering the 
convolution of different exponential waiting time distributions. In this case we have 
f — fi * fii where each fi is of the form Eq. (5) with parameter Aj, which leads to 

f(0 = 2- ■ ; -=- — : ; 7 < 1 4) (23) 



s v /rTr 4f 2 V s 



where s, p denote sum and product of the parameters Aj. Setting Ai = A and A 2 = rX, 
so that the parameter r = A2/A1 gives the relative scale among the rates of the two 
waiting time distributions the function q(t) takes the 



q(t) =e 2 



cosh f \/ r 2 — 6r + 1 — ) H = = sinh ( \Jr 2 — 6r + 1 — 

2 J Vr 2 - 6r + 1 \ 2 



,(24) 



which can oscillate and take on negative values, thus leading to a non zero measure of 
non-Markovianity, for 3 — 2y / 2 ^ r ^ 1/(3 — 2y / 2)- Also in this case non-Markovianity 
does not appear when one of the rates is much stronger than the other, so that the overall 
distribution is dominated by only one of the two components, which is memoryless 
distributed. The sign of q(t) as a function of the relative strength r and the rescaled 
time At is plotted in Fig. 2. Note that also in this case the measure of non-Markovianity 
is given by the sum of the maxima of \q(t)\. 
The value A = (0, |, \, 0) leads to the map 

1 1 

£ P P= TfxPVx + 2 a yP a y = a +P a - + °"-P°"+> 

for which also populations are affected, since fi = (1,0,0,-1) and A = 
(1, g (t) , g (t) , q (t)), with g (t) survival probability as in Eq. (7). As discussed in [9] 
the measure of non-Markovianity of this map according to the trace distance criterion 
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Figure 2. (Color online) (left) Plot of the sign of q(t) given by Eq. (24) for the 
convolution of two exponential distributions with Ai = A and A 2 = r A as a function of 
r and Xt. The black regions correspond to a negative sign, (right) Plot of the sign of 
Ai-2i/(t) for the same waiting time distribution as a function of v and At. One clearly 
sees the threshold for the inset of non-Markovianity at v = 1/2. 



is the same as for the dephasing map characterized by A = (0, 0, 0, 1), so that in both 
cases it is determined by the waiting time distribution. As we shall see this is no more 
true for the other measure of quantum non-Markovianity. 

In the case A = (1 — u, 0, 0, is), with ^ v ^ 1, one has a mixture of 
phase-flip and identity, leading to fi — (1, 1 — 2v, 1 — 2u, 1) and therefore A (£) = 
(1, \i-2v (t) , \i-2v (t) , 1), with \i-2v (u) as in Eq. (19). The resulting map again 
describes pure dephasing, but the off-diagonal matrix element is multiplied by a factor 
Ai-2i/ (t). In order to assess non-Markovianity one now has to study \i-2u (t) instead of 
q (t). Considering again the same class of waiting time distributions, one sees that now 
a threshold appears, so that one can have non-Markovianity only if v > |, as shown 
in Fig. 2. This shows that indeed non-Markovianity depends on both elements of the 
couple {£, f (t)}, at variance with the classical case and the previous examples. 

The criterion of non-Markovianity introduced in [8] relies on the requirement of 
divisibility of the quantum map in terms of completely positive maps, so that Eq. (3) 
applies where now each A is a completely positive map. To detect the failure of 
divisibility one can consider the Choi matrix associated to the map A (t, s), with t ^ s, 
which can be obtained from Eq. (20). For the general case the expression is somewhat 
cumbersome, and an alternative way is to look at the sign of the coefficients in the 
master equation in time-convolutionless form corresponding to Eq. (13). For the case 
A = (0, |, |, 0) this path has been followed in [9], the master equation reads 

dp 1 ( \ x (t) Xy (t) Mt)V lA*(f) ^ ! 

with X z (t) = q(t), X x (t) = X y (t) = g (t). In this setting one can consider situations 
in which the measure of non-Markovianity related to the trace distance is zero, since 
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the map is divisible in the sense of positive maps, but the intermediate maps are not 
completely positive. This is the case for the convolution of two different exponential 
waiting time distributions as in Eq. (23), for a suitable choice of the ratio A1/A2. Note 
that for the situation in which all A« (t) are different one has 

dp 1 ( K (t) X y (t) X z (t)\ lX z (t) 




A, it) 


,\ (t) 


X x (t) 


\ it) 


X x (t) 


Xy (t) 


X x it) 


Xy (t) 



1 X X (t) Xy(t)\ 

[a * pa * ~ p] ~i [xaT) ~ \W) ) [(Typ(Ty " p] ' 

and the two intermediate channels have opposite coefficients, so that unless X x (t) = 
X y (t), which is the case considered in [9], one is always negative. In this case however 
one cannot read divisibility from the sign of the coefficients, since the Lindblad operators 
appearing in them are not linearly independent, indeed Eq. (25) can be written as 

dp ifj,W K(t) Mt)V„„ 1 (9fi > 

s= + 4^M^)-M^)-M^)J l ' 7 ^- ^ ' , (26) 

l(K(t) \(t) K(f)\. 1 (\ x (t) X 9 (t) k,(t)\, . 

In the present framework we can indeed point to a situation in which, due to a subtle 
balance, all coefficients in Eq. (26) are positive, thus granting divisibility, but this is no 
more true for the coefficients of Eq. (25). This situation is depicted in Fig. 3 considering 
A = (|,|,|,|) and a suitable waiting time distribution. Even though this fact can be 
easily understood from a conceptual point of view, it is useful to stress it by means of 
an explicit example. 

For the case A = (0, 0, 0, 1) of pure dephasing in order to detect non-Mar kovianity 
it is convenient to consider the Choi matrices associated to the maps A(t + s,t) = 
A(t + s, 0)A _1 (t, 0), so that according to Eq. (20) we obtain 

A> kl (t + s,t) = Mag (l + «&±^, 0,0,1- «£±^) , (27) 

where the simplicity of the result strongly depends on the convenient choice of basis. 
Failure of divisibility is then detected when at least one of the coefficients of these 
collection of matrices depending on two temporal indexes becomes negative. The sign 
of the smallest non zero eigenvalue is plotted in Fig. 3 for a waiting time given by 
the convolution of two identical exponential distributions. In accordance to the result 
obtained relying on the criterion based on trace distance distinguishability, the map is 
indeed non-Markovian. Note however that the divisibility property does depend on the 
initial time considered, i.e. for certain time windows, which obviously include the initial 
time t — 0, the maps A(t + s,t) are completely positive for any time interval s. At the 
same time the violation of complete positivity does decrease for large s. Moreover while 
the two criteria agree in labelling the map as non-Markovian, as discussed in detail in [9] 
they assign to it different measures. In particular as follows from Eq. (27) the approach 
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Figure 3. (Color online) (top left) Plot of the two distinct coefficients of Eq. (26). 
(bottom left) Plot of the four distinct coefficients of Eq. (25). In the latter case 
one is always negative. The coefficients of the two equivalent expressions of the 
time-convolutionless master equation associated to the map Eq. (21) for the case 
A = (|)§,|,|), so that (j, = (l,|,-i,-|) are obtained from Eq. (17). The 
considered waiting time distribution is of the form Eq. (23) with A2 = r\ and 
r = 0.13. (right) Plot of the sign of the lowest non zero eigenvalue of the Choi 
matrix corresponding to the map A(t + s, t), for a stochastic map corresponding to 
pure dephasing and waiting given by the convolution of two equal exponentials. The 
black regions correspond to a negative sign. For certain values of t, corresponding to 
the points in which the inverse does not exist, in the plot denoted by a gray vertical 
line, the memory is arbitrary long, while it stays finite for other values. 



based on divisibility assigns an infinite measure to this map as soon as the quantity q(t) 
goes through zero, so that at variance with Eq. (22) processes with different memories 
are put on the same footing. Let us note that this difficulty can be overcome by keeping 
divisibility of the time evolution in terms of completely positive maps as a signature 
of Markovianity, quantifying however its violation in a different way As suggested in 
[34] one can introduce a different weight, considering the integral of the arcotangent of 
the sum of the negative eigenvalues of the matrix Eq. (27) in the time regions in which 
complete positivity breaks down, renormalizing by the extension of these regions. For 
the case at hand indeed this modification makes the non-Markovianity measure based 
on divisibility finite, even though the expression Eq. (22) for the measure based on 
distinguishability remains much easier to evaluate. 



4. Conclusions 



In the present manuscript we have considered a class of quantum dynamics which can 
be obtained from a generic unital stochastic map on C 2 and a classical waiting time 
distribution, thus extending the class of examples considered in Ref. [9] to show how 
the recently introduced notions of quantum non-Markovianity relate to the classical 
one. The considered examples show how versatile the class of semi-Markov processes 
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and their quantum counterpart can be in order to study the notion of non-Markovian 
process in the quantum framework and highlight different possible behavior. 
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